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Abstract 

We study dynamics of flat directions in the minimal supersymmetric standard model taking 
account of its constituent fields. It is found that there exist new instabilities due to the D-term 
potential and the nature of these instabilities depends on the eccentricity of the orbit. For a 
roughly circular orbit, it is similar to narrow-band resonance described by the Mathieu equation. 
For an elliptic orbit, the instabilities exhibit tachyonic nature. In the narrow-band regime, we show 
that the existence of the instabilities is crucial to the formation of Q-balls. The Q-ball formation 
proceeds through quasi-stable state called "/-ball" . The transition from /-balls to Q-balls is made 
efficient by the D-term instability and during this process some fraction of the charge of the Q-ball 
is emitted. This discovery may revive the scenario which relates the baryon number and dark 
matter density of the universe. Furthermore, the tachyonic Z?-term instability causes the drastic 
decay of the flat direction moving in an orbit with relatively large eccentricity. Thus the evolution 
of flat directions is completely altered by the appearance of this instability. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



The scalar potential in the minimal supersymmetric standard model (MSSM) has flat 
directions along which there is no classical potential in supersymmetric (SUSY) limit. In 
the real world, the potential along the flat directions is not completely flat because SUSY 
breaking induces soft masses of the order ~ 0(100) GeV. Such flat directions play important 
roles in cosmology. In the inflationary universe some fields corresponding to flat directions 
have large field values during inflation, and start oscillation after their mass becomes com- 
parable to the Hubble parameter. The flat directions generally consist of squarks, sleptons 
and Higgs, and hence have baryon and/or lepton numbers. Thus, the oscillation of the flat 
directions and subsequent decay can produce large baryon and/or lepton numbers in the 
universe if some baryon (lepton) number violating terms exist in the potential. This mech- 
anism was discovered by Affleck and Dine 0] and it has been one of promising scenarios 
for explaining the baryon asymmetry of the present universe. Hereafter we call scalar fields 
corresponding to flat directions Affleck-Dine (AD) fields. 

Recently, it was found that condensates of some AD fields evolve into non-topological 
solitons called Q-balls through strong instabilities during oscillation 0, 0|. In fact, nu- 
merical simulations showed that initial tiny fluctuations of the AD field rapidly grow and 
deform into large Q-balls jH[ 0| . The existence of Q-balls complicates dynamics of the AD 
field and significantly changes the scenario of the AD baryogenesis jjjl 0, 0, H| . Q-balls also 
provide several interesting possibilities. For example, Q-balls can account for the dark mat- 
ter and baryon asymmetry of the universe simultaneously 0, H| , or they may produce very 
large lepton asymmetry jj| E| • 

So far, dynamics of flat directions has been studied with use of "single-field parametriza- 
tion", i.e. a flat direction is described by one complex scalar field. In fact, however, the flat 
direction is actually composed of several scalar fields. When one only considers the dynamics 
of a homogeneous AD field, the single-field parametrization may be good enough. However, 
we now know that the fluctuations of the AD field is crucial to Q-ball formation. Since the 
field fluctuations increase the number of degrees of freedom in the dynamical system, the 
single-field parametrization is not always adequate to describe the evolution of the fluctua- 
tions and multi- field analysis is necessary. In the previous paper we demonstrated the 
importance of the multi-field analysis when one of component fields that form Q-balls has a 
large decay rate. 

In this paper, we study the dynamics of the AD field taking account of its constituent 
fields. We find that there exist instabilities due to the D-term potential and the nature 
of these instabilities is determined by the eccentricity of the orbit. For a roughly circular 
orbit, it is described by narrow-band resonance. However, as the eccentricity increases the 
instability exhibits another character: tachyonic instability. 

In the former case, since the instability is the narrow-band resonance, its effect on the 
cosmo logical evolution of the homogeneous mode of the AD field is not so important. How- 
ever, the existence of the instabilities is crucial to the formation of Q-balls. In general, in 
gravity mediated SUSY breaking models, the Q-ball formation is complete through quasi- 
stable state "/-balls" 0] inside which the AD field rotates in an elliptical orbit. It is found 
that the D-term instability quickly makes the orbit of the AD field circular and a Q-ball con- 
figuration is reached. During this process some fraction of the charge confined in the Q-ball 
is emitted, which is different from the previous result obtained with use of the single-field 
parametrization. As a result the baryon density of the universe is accounted for by both 
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Q-ball and residual baryon emitted during the Q-ball formation. On the other hand the 
dark matter is explained by the lightest SUSY particle (LSP) produced in the Q-ball decay. 
This may revive the scenario which relates the baryon number and dark matter density of 
the universe. 

Meanwhile, if the instability becomes tachyonic, its effect on the cosmological evolution 
of the AD field is no longer negligible. Irrespective of whether Q-balls are formed or not, 
the homogenous motion always suffers from rapid grow of the instabilities. Therefore the 
AD field moving in an orbit with large eccentricity decays soon after it begins to oscillate. 
We expect that the similar instabilities will exist in many other cases. For instance, the 
instabilities might change the reheating processes of some supersymmetric inflation models. 

In the next section, we investigate the D-term instability, and the effect on Q-ball for- 
mation is studied in Sec. IIHI Sec. |IV] is devoted to the discussion on the fate of the flat 
directions under the effect of the tachyonic instability and finally we give discussions and 
conclusions in Sec. |V| 



II. D-TERM INSTABILITY 

In this section, first we explain how flat directions are expressed, paying particular atten- 
tion to the relationship between single-field and multi-field parametrizations. The former 
is usually adopted in many literature, since it well describes the collective motion of the 
constituent fields. Interestingly however, for instance, if some of the constituent fields decay 



into lighter species, such simple description fails 11]. Furthermore, the D-term instability 
appears only in the multi-field parametrization, which enables us to explicitly deal with the 
D-term potential. After reviewing the flat direction, we will investigate the instability in 
detail. 



A. Flat directions 

In the MSSM, there are many flat directions along which both the D-term and F-term 
potentials vanish at the classical level. The D-flat direction is labeled by a holomorphic 
gauge-invariant monomial, X. Most of the D-flat directions can be also F-flat, simply due 
to a generation structure of the quark and lepton sector. Since it is easy to satisfy the F-flat 
condition, let us concentrate on the F-flat condition in the following. The F-flat direction, 
X, can be expressed as 

N 



x=n$i, (i) 



where N superfields constitute the flat direction X, and we have suppressed the gauge 
and family indices with understanding that the latin letter i contains all the information to 
label those constituents. When X has a nonzero expectation value, each constituent field 
also takes a nonzero expectation value 

<*<> = ^f 6i - (2) 
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Here each <pi/y/2 is the absolute value of the expectation value and is related to every 
other due to the D-flat conditions: 



where T4 are hermitian matrices representing the generators of the gauge algebra, and 
they are labeled with A. This condition can be usually satisfied if we take all (pi equal: 
(pi = <j)/y~N. Note that the D-flat condition dictates that all the amplitudes be equal, 
while the phases, 9i, remain to be arbitrary. In fact, it is necessary to know the pattern 
of baryon and lepton symmetry breaking, in order to identify the Nambu-Goldstone (NG) 
boson relevant for the AD mechanism. The situation becomes simple if the spontaneously 
broken symmetry coincides with the explicitly violated one. This is the case if the A-term 
is given by some powers of X. Then all 9i become equivalent, and the NG boson can be 
identified with the average of 9i j3|: 

i 

where 9 is orthogonal to NG modes corresponding to the spontaneously broken £7(1) gauge 
symmetries: the linear combination of left isospin Tlz and weak hypercharge Y . This is 
because ^ Tl3,i = ^ Y% = by definition of the Z)-flat direction. Thus, the dynamics of 
the flat direction, X, is described by one complex scalar field, $: 

* - fc. <*) 

Note that $ is canonically normalized. The single-field parametrization is so useful in usual 
situations that we almost forget that the flat direction is actually composed of multiple 
fields. However, one must keep in mind that this parametrization is just an approximation 
which becomes exact only in the limit of D-flatness, and that it describes only the collective 
motion of the constituent fields. The smallness of the D-term potential does not necessarily 
guarantee that its effect on the dynamics of flat direction as well is negligibly small. This is 
actually what we are going to show in this paper. 

Flat directions are lifted by the SUSY breaking effect and non-renormalizable terms. To 
be concrete, let us adopt the gravity-mediated SUSY breaking model. The flat direction is 
then lifted as 

v s (m) = E^n 2 ( 1+ ^ 1o s(^)) ( 6 ) 

where m ; is a soft mass, Ki a coefficient of the one-loop correction, M* the renormaliza- 
tion scale to define the mass. This potential reduces to the following in the single-field 
representation, 

V.(*) = m!|$| 2 (l + inog(|^), (7) 
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where m$ and K are denned as 

— Vm 2 K — ^i m * Ki (g) 
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The absolute value of K is estimated as \K\ = 0.01 ~ 0.1 0. Moreover, assuming a 
nonrenormalizable operator in the superpotential of the form 

X k 

Wnr ~ kM Nk ~ 3 ' ^ 
the flat direction is further lifted by the potential 

N 



vMm) = Ei%p> ( 10 ) 

i=l 

where M is a cutoff scale. In fact, the nonrenormalizable superpotential not only lifts the 
potential but also gives the baryon and/or lepton number violating A-term of the form 

V *({*<}) = "rnj^^X* +h.C., (11) 

where m 3 / 2 is the gravitino mass, a m is a complex constant of order unity, and we assume a 
vanishing cosmological constant. 

During inflation, the flat direction, X, is assumed to take a large expectation value. After 
inflation ends, the Hubble parameter starts to decrease, and becomes comparable to the mass 
of $ at some point. Then the flat direction starts to oscillate and acquires a finite angular 
momentum due to the A-term. Since we are interested in the evolution of the flat direction 
after it starts to oscillate, we neglect non-renormalizable terms (including the A-term), and 
take ellipticity of the trajectory as a free parameter. Instead of eccentricity, however, we 
use axial ratio e defined as the ratio of the magnitudes of the major axis and minor axis of 
ellipse described by the AD field (see Eq. (J53J)). For instance, e = 1 corresponds to a circular 
orbit, while e = represents a straight-line motion. In the following sections we study the 
evolution of with the potential Vc>({$i}) + K({^i}) m the multi-field parametrization, 
while the potential is simplified as Eq. (JJJ) in the single-field parametrization. 



B. Instability associated with D-term potential 

Now let us consider dynamics of a flat direction taking account of the D-term potential. 
As a beginning we would like to clarify a setup relevant for our discussion. We take up the 
simplest possible flat direction composed of two scalar fields with U(l) gauge symmetry, 
X = $i(+l) $ 2 (— 1), but it is trivial to extend our results to more generic case. The effect 
of the expanding universe is neglected for the moment, but we will get back to this point 
later. 

The constituent fields, $1,2, obey the equations of motion 

, ( + ^M. , (j = li2) (12) 
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with 




where the overdot denotes the differentiation with respect to time, g is the U(l) gauge 
coupling constant, and we take K\ = K2 = to concentrate on the instability that originates 
from the D-term. Rewriting the equations of motion in the real and imaginary components: 



we obtain 




for the real parts. The equations for the imaginary parts can be obtained by interchanging 
the subindex R with I. When the flat direction starts to oscillate, its amplitude is usually 
much larger than m$/g, so that the D-flat condition: |$i| = |$ 2 |, is satisfied with good 
precision. With an appropriate phase rotation, the D-flat condition can be also expressed 
as 4>ir = 4>2R and 4>u = 4>2i- Thus the homogeneous motion can be described by a single 
field $ rotating in a parabolic potential. 

Let us now consider fluctuations of the flat direction moving in an elliptical orbit with 
axial ratio 1 : < e < 1. First we derive the equations of motion for linearized fluctuations 
from Eqs. (|15j) and Using the .D-flat condition of the homogeneous motion, we have 

54>ir + (k 2 + m\) 54>ir + g 2 4>iR {(j)x R 54>m + (puScpu - 4>2rS4>2r - faiSfoi) = 0, (17) 

54>2R + {k 2 + mj) 5<p2R + g 2 4>2R {4>2R<><j)2R + <p2I^<p2I ~ 4>IR$<I>1R ~ <PllHll) = 0, (18) 

and similar equations for the imaginary parts. Here and in what follows we consider fluctu- 
ations in momentum space, and k denotes the momentum. Further simplification is possible 
if mi = m 2 {= m$), so we assume that this is the case just for simplicity. However, note 
that this assumption is not obligatory, and that it does not change generic features of the 
D-term instability. From Eqs. (fT7j) and (|TH|) . one can see that 54>ir + 54>2R obeys a simple 
harmonics equation with a constant frequency. In other words, 5<p\R + 5<p2R does not grow 
at all. Similar arguments apply to the imaginary parts 54>u + 5<p2i as well. Therefore, for a 
growing mode, it is a fairly good approximation to take 

HlR + S(j)2R = 0, 5<t>u + 5(j)2I = 0. (19) 



We do not consider the case that e is extremely close to 0, which corresponds to an almost straight-line 
motion. This is because other type of instability known as "preheating" |l4| comes into play. In particular, 
the instant preheating could proceed if e < m$ / g|$| osc , where |$| OS c represents the amplitude of 

the flat direction when it begins to oscillate. Since the flat direction generally acquires a finite angular 
momentum due to the A-term, we can easily avoid such drastic decay. 
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Note that this approximation picks out the relative fluctuations between the constituent 
fields $1 and $2, which are disregarded in the single-field parametrization. See Fig. ^ If all 
the relative fluctuations possessed very heavy masses ~ g ($), the single-field parametriza- 
tion might suffice. However, one of the relative fluctuations is actually a massless mode as 
shown below, which experiences the D-term instability and thereby grows rapidly. 
Hereafter we follow the motion of relative fluctuations 5(f>R and 5(f>i defined as 



Hr - 5(f> 2 R 

2 

hi ~ H21 



(20) 



and drop the subindices "1" and "2" at the homogeneous fields. The equations of motion 
for fluctuations 5<f) R and 5(f>i are written in the following form: 



J ^ 3?'=° (21) 



with 



2 _ / k 2 + m% + 2g 2 4> 2 R 2g 2 <p R <p I , 
Uk ~ I 2g 2 (f> l( f> R k 2 + m% + 2g 2 (f) 2 > " 



In order to diagonalize we define and 5(f>M,m as 

■h)- Rk \^ m )' Rk = J \h <Pr J' (23) 

where <p is a constant representing the typical amplitude of the flat direction. The eigen- 
values of io\ corresponding to 5(f) m and 5(f> m are k 2 + m\ + 2g 2 ((p 2 R + (ft 2 ) and k 2 + m|, 
respectively. Since the frequency of 5(f) m is much larger than that of the homogeneous mo- 
tion, its evolution is adiabatic, which means that the particle production of 5(f>M is negligible. 
The equation of motion for 5<f) m is thus given by 

sin 2?n<f>£ 

5(f) m - m*(l - e 2 ) — — - 2 -5(f) m + k 2 5(f) m = 0, (24) 

cos z m$t + e z sin m$t 

where we have substituted the solution for the homogeneous mode: 

4>r — 0o cos m$i, 

(pi = E(f> sin m$t, (25) 

and we have neglected 5<pM in Eq. (J23|) since it is rapidly oscillating. 50 m is fluctuation 
along the elliptical orbit as shown in Fig. |21 which can be confirmed by noting that Eq. (J2%Jl 
except for the third term is same as the equation of motion for the angular variable 9 of the 
homogeneous AD field moving in an elliptic orbit with axial ratio e: 



6 + 2-j6 = 0, (26) 



where (f> and 9 represent the amplitude and angle of $ in the single-field representaion: 
4> 2 /2 = <f) 2 R + (f> 2 ; tan# = 4>i/4>r- Using this analogy, it is easy to see that Eq. (J23j) contains 
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instabilities, leading to the rapid growth of 5<j) m . Clearly, the angle 9 increases monotonically 
as the AD field moves around in its orbit. In particular, 9 exhibits step- function-like behavior 
for £ <C 1. Therefore we can deduce that 5<p m as well increases monotonically aside from 
the case of e = 1, in which 5<j) m oscillates with the constant angular frequency k, while 9 
increases linearly with time. 

To see the band structure of this instability, let us rewrite this equation in terms of r and 
5ip defined as 

r = m 9 t- |, (27) 



dip ee 6<t> m y/l + e 2 - (1 - e 2 ) cos2r. (28) 
Then we obtain the final result: 

5<p" + ( + F(e, t)) dp = (29) 



with 



_ (1 -6- 2 )(sin 4 r-e 2 cos 4 r) 

F{e,r) = , — , (30) 

sin r + cos z rr 



where the prime denotes the differentiation with respect to r. It is worth noting that the 
amplitude of the flat direction <p disappears in the Eqs. (J2*%jl and (|29|). Therefore the 
instability band and exponential growth rate are independent of the amplitude. 

First let us take the limit of e — > 1 in Eq. (|29|). It then reduces to the well-known Mathieu 
equation fvi\ : 

dp" + (A - 2q cos 2r)5p = (31) 

with A = k 2 /m% and vanishingly small g-parameter q ~ 1 — e. For q <C 1, the resonance is 
known to occur only in narrow bands near A = I 2 , (I — 1, 2, . . .), and the instability vanishes 
if q = 0. Therefore we expect that the nature of the instabilities of Eq. (J2U|) is similar to that 
of the Mathieu equation, if the AD field is rotating in a roughly circular orbit. That is to 
say, the narrow-band resonance occurs around k/m<$, — I, (I — 1, 2, . . .), and 5<p grows unless 
the flat direction is moving in an exact circular orbit (i.e., q = 1 — e = 0). The instability 
chart of Eq. is illustrated in Fig. El and one can see that the narrow-band structure 
around e — 1 resembles those of the Mathieu equation. 

Next we take up the other limit of e 1. In this case, the function F(e,r) behaves as 
follows. For \t\ < £ -1 / 2 , it becomes negative and the minimum locates at r = 0: F(e,0) = 
—e~ 2 , while F(e,t) — 1 for |r| > e~ l l 2 . See Fig. 0] Therefore those fluctuations with 
k < e -1 m$ experience tachyonic instabilities and grow exponentially each time r passes 
t = nrr, (n = 0, 1, 2, . . .). That is to say, the instability shows the step-function-like behavior, 
as expected. As can be seen from Fig. 01 the instabilities become stronger for smaller e. 
However we expect that the exponential growth rate will converge to some finite value 
~ 0(1) as e -> 0. 

So far we have neglected the cosmic expansion, but the above-mentioned nature of the 
instability does not change much even if the cosmic expansion is taken into consideration. 
For simplicity, let us assume that the universe is dominated by the energy of the oscillating 
inflaton, when the AD field starts to oscillate (i.e., H ~ m$). Then similar arguments lead 
to the equation of motion of S(f) m given by Eq. (J24|) with k 2 replaced with k 2 /a 2 , where a 
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denotes the scale factor. Note that the solution of the homogeneous mode is now given by 



4>R = 4>oQ> 3 ^ 2 cosm$t, 

<pi = e<f) a~ 3 / 2 sin m$i. (32) 

Due to the redshift of the momentum k, the band structure becomes vague as seen in Fig. |SJ 
However, except for e ~ 1 where the instability is narrow-band resonance, there still exist 
strong D-term instabilities because of their tachyonic nature. 

The above argument is not limited to the parabolic potential, but it can easily be gen- 
eralized to an arbitrary potential with U(l) symmetry. Although the collective fluctuations 
as well are subject to instabilities peculiar to the potential in general, we can concentrate 
on the relative fluctuations by adopting Eq. (fTUjl since they are independent. As long as the 
homogeneous motion satisfies the D-flat condition, the lighter one of the relative fluctuations 
satisfies 

50 m + 2^50 m + k 2 5<p m = 0, (33) 

where the definitions of the variables are same as before. If we rewrite this equation in terms 
of 5<p = <50 m , we have 

5(p + ^k 2 - | j 5tp = 0. (34) 

When the AD field comes the closest to the origin, 0/0 is expected to take its maximum 
value, which leads to the tachyonic instabilities for k 2 < [0/0] max- Also if the homogeneous 
motion in the limit of e — > 1 can be approximated as 

<p R = 4> Q cos Qt, 

(pi = e(j) smflt, (35) 

then Eq. (J34j) reduces to the Mathieu equation with A = k 2 and q ~ Q 2 (l — e). In this case, 
the instability bands are narrow and locate at k/Q — I, (1 = 1,2,.. .). Thus the existence 
of the D-term instability does not depend on the specific form of the potential. 

In this section we have considered a simple model of D-flat directions. However the 
existence of the instability should be common to all D-flat directions. The reason is as 
follows. First, all the MSSM flat direction is at least charged under U(l)y- Thus D-term 
potential similar to Eq. (|13j) always exists. Second, the instability originates from a massless 
mode of the relative fluctuations between the constituent fields, the existence of which is a 
common feature of D-flat directions. To sum up the major characteristics of the D-term 
instability, (i) for e ~ 1, it is narrow-band resonance similar to the one expresssed by the 
Mathieu equation; (ii) the tachyonic instability appears for e < 1; (iii) it is independent 
of amplitude of the flat direction. In the following sections we will see how such D-term 
instability affects the dynamics of the AD field. 



III. EFFECTS ON Q-BALL FORMATION 
A. Review of Q-balls and /-balls 

Let us now discuss how the D-term instability affects the Q-ball formation processes. 
Before we go into details, it will be useful to briefly review the Q-ball formation. For the 
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FIG. 1: The relative fluctuations between $1 and $2 (solid arrows), which grows due to the D- 
term instability. The single-field parametrization deals only with the collective fluctuations shown 
as dotted arrows. 
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FIG. 2: The homogeneous motion and the fluctuation along the elliptical orbit. 




FIG. 3: The instability charts of Eq. (|29[). The momentum fc is normalized by m$. The left and 
right panels are equivalent except for the plotted range of e. Each contour line corresponds to 
fjb/m^ = 0.1,0.3,0.5,0.7,0.9 from right to left in the left panel, while /x/m$ = 0.1,1,2 from right 
to left in the right panel, where [i is an exponential growth rate of fluctuations. 
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FIG. 5: Same as Fig.|3]except that the cosmic expansion is taken into account. We take the initial 
condition as H = m$, and assumed the matted-dominated like universe. 



moment we adopt the single-field parametrization. Let a flat direction $ possess an unit 
baryon charge. During inflation $ is assumed to take a large expectation value. When 
the Hubble parameter becomes comparable to the mass of $ after inflation ends, $ starts 
to oscillate in the potential Eq. (|7jl. If the sign of K is negative, it experiences spatial 
instabilities, leading to the Q-ball formation U 0, 0|. The profile of thus formed gravity- 
mediation type Q-balls is approximated by the spherically symmetric gaussian form with 
good precision: |4( 



Mr,t) 



>r + i4>i 
V2 ' 



e 
e 



~ r2 / R Q cos cut, 
- v2 / r q sin cut, 



(36) 
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with 



uo = m^\/l + 2\K\, 
R Q = -7^—^ (37) 



where and Rq denote the angular velocity and radius of the Q-ball, respectively. 

Although a Q-ball solution is realized as the energy minimum state with a fixed global 
U(l) charge, it is a /-ball that is first formed in general. The profile of /-balls is almost 
same as that of Q-balls except for axial ratio of the orbit: 

<j> R (r,t) = <p e' r2/R Q cos tut, 

7 (r,t) = e I( p e- r2/R Q smut, (0 < ej < 1). (38) 

The reason why /-balls are formed instead of Q-balls is that the quadratic potential Eq. (J7J) 
allows the dynamics of system to have another invariant, the adiabatic charge / and 



that the axial ratio of the homogeneous motion is generally less than 1. Thanks to the 
conservation of the adiabatic charge, /-balls are formed as a quasi-stable state, and the 
axial ratio of the orbit inside /-balls, Ej, is almost same as that of homogeneous mode before 
fluctuations grow, Eh- 

El ~ E H . 

However this relation does not hold when the D-term instability is taken into account, and 
Sj is larger than or equal to Eh as shown below. 

Central to this issue is how /-balls eventually transform into Q-balls. Due to the elliptical 
orbit, /-balls are considered to be an excited state of Q-balls. However it was unclear to 
date how /-balls emit the extra energy to settle down to Q-balls, since they have very long 
lifetime, which makes it difficult to follow the evolution in numerical calculations. Here we 
would like to argue that the D-term instability violates the conservation of the adiabatic 
charge, which enables /-balls, if formed, to immediately transform into Q-balls by emitting 
both the energy and U(l) charge. 

Before we proceed, however, it is necessary to ask whether /-balls are really formed 
under the influence of the D-term instabilities. Since the exponential growth rate of the 
/-ball formation is given by fig ~ |/T|m$ ~ 0.1m$ [4], the /-ball formation will complete 
before the D-term instabilities grow sufficiently, if the initial axial ratio Eh is larger than 
~ 0.7 (see Fig. |SJ). Furthermore, even if the D-term instabilities grows faster than ~ |/f |m, 
some /-ball-like objects might be generated as a result of the interactions between the 
developed fluctuations, because the fastest growing mode of the D-term instability locates 
around k/a ~ m$ which is not so different from the typical size of /-balls. In fact, we have 
performed numerical calculations not taking account of the cosmic expansion and confirmed 
that /-balls are formed for eh > 0.1, but the axial ratio of the orbit inside thus formed lumps 
is larger than eh from the beginning, e.g., ej ~ 0.3 for eh = 0.1. In particular, /-balls with 
Ej < 0.2 are not formed, probably because the D-term instability with Ej < 0.2 is too strong 
for /-balls to retain their configuration. Although the cosmic expansion might make the 
interactions between the fluctuations less effective, some part of the charge asymmetry will 
still be stored inside such /-ball-like objects. We will discuss the evolution of the axial ratio 
of the trajectory inside these objects in the next two subsections. 

The /-ball formation is not the end of the story. Once /-balls are formed, the dynamics of 
the AD field inside them decouples from the cosmic expansion. Since the D-term instabilities 
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lie at smaller scales than a typical size of /-balls, fluctuations inside them continue to grow 
exponentially. It should be noted that the D-term instability does not depend on the 
amplitude, so that they will grow throughout the /-ball. The growth will last either until 
the conversion of /-balls into Q-balls completes (i.e., the axial ratio e becomes 1), or until /- 
balls are completely destroyed by the instability. We need to resort to numerical calculations 
to determine which is realized. 

B. Numerical Calculations 

As mentioned in the previous subsection, we have numerically investigated the formation 
processes of /-balls taking account of the D-term instability and found that almost all the 
charge asymmetry is taken in the /-balls for Eh, in\ > 0.1: 

// * 1, (39) 

where // denotes the fraction of the charge absorbed by the /-balls. Due to the technical 
reasons, it is difficult both to take account of the cosmic expansion and to follow the evolution 
of the AD field with SH,\m < 0.1. The cosmic expansion might change the /-ball formation 
inefficient to some extent. The strong instability for eujm 0.1 as well prevents the /-balls 
to be formed. But what we would like to show here is that the charge asymmetry is partly 
left behind during the Q-ball formation, even if it is once absorbed by /-balls. 

Assuming that /-balls are formed, we have followed the evolution of the /-ball in order to 
show that the D-term instability actually occurs inside /-balls. We have solved the equations 
of motion given by Eq. (|12|) on the one dimensional lattices with the initial condition: 

4>ir = 4>2R = (j) e~ r2/R Q cosm$t, 

01/ = 02/ = £/ O e" r2/fl « sinm$t. (40) 

Once /-balls are formed, the dynamics inside them decouples from the cosmic expansion. 
Therefore we can safely neglect the effect in contrary to the case before the /-balls are formed. 
It should be reminded that the axial ratio £/ in the above equation does not necessarily 
coincides with that of the homogeneous motion eh- The space x and time t are normalized 
by the mass scale m$, while the field value is normalized by M*. We take the following 
values for the model parameters: vn\ — mi — m<j>, 4> /M* = 1, g 2 Ml = 10 3 or 10 4 m|, and 
K\ = K 2 = —0.1. The typical evolution of the axial ratio £j with the initial value: Ej = 0.5, 
is illustrated in Fig. |UJ The spatial distributions of the energy and the charge are shown in 
Fig. [7| As seen in the figure, the axial ratio increases, but it does not reach 1 due to the 
following reasons. First of all, as the axial ratio comes closer to 1, the instability becomes 
weaker as shown in Fig. El Thus in the finite CPU time the complete conversion of the /-ball 
to a Q-ball cannot be seen. Second, the finite lattice size allows the decay product emitted 
from the /-ball to be absorbed by the /-ball again. This effect as well prevents the /-ball to 
transform into the Q-ball. (We have chosen the time interval in Fig. so that this effect is 
negligible.) However, in the real world, such a tendency toward Q-balls will last throughout, 
because the adiabatic charge is no longer conserved due to the D-term instability. We thus 
expect the conversion will be completed relatively soon after the /-ball formation. 

We have also calculated the power spectra of one of the constituent fields, $i. See Fig. |SJ 
We take the size of the /-ball larger (i.e., smaller \K\) to obtain clear signature of the 
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instability band. The result coincides with the instability band obtained by the analysis of 
linearized fluctuations. Thus we conclude that the D-term instability plays an essential role 
to transform /-balls into Q-balls. 



C. Modified relation between baryon and dark matter abundance 

What should be emphasized at this juncture is that both the energy and charge are 
emitted from the /-ball. This is because the four-point interaction in the D-term inevitably 
induces both the charge-transferring (e.g., — ► $i$i) and charge-conserving processes 
(e.g., — > $!$*). If almost all the charge is absorbed by the /-balls at their formation, 
the fraction of the baryon charge stored inside Q-balls, /q, depends only on the initial axial 
ratio 2 . Since /q is defined as the fraction of the charge stored inside Q-balls assuming all 
the charge asymmetry is once absorbed by /-balls, the actual fraction kept inside Q-balls 
is fs = fifq- It is important to know the value of fs to relate the baryon number to 
dark matter density. The numerical results of fg for several values of Ej are summarized 
in Table |U The values of Jq shown in the table should be considered to give upper bounds 
of Jbi since some part of the charge is left behind at the /-ball formation and also the 
conversion is not yet completed. It should be noted that the /-ball comes apart for Ei < 0.2 
due to the strong D-term instability. This fact agrees with the result that /-balls with the 
axial ratio smaller than ~ 0.2 are not formed. 

The neutralino LSPs produced by the decay of Q-balls survive and can be the dark 
matter, if the decay temperature of Q-balls is lower than the freeze-out temperature of the 
LSP. Let iV be the average number of the LSPs produced by the decay of unit baryon charge 
of Q-balls. The abundance of the LSP dark matter is then related to that of baryon number 
as follows. 



where £Idm,b denote the density parameters of the dark matter and the baryon, m L sp and 
m n are masses of the LSP and the nucleon, respectively. Substituting the WMAP results |18( : 
VL DM h 2 = 0.116 and Vt B h 2 = 0.024, we obtain 



Previously, /s ~ 1 was supported by the numerical calculations [6j, therefore the coincidence 
of the baryon and dark matter abundance cannot be explained in this way. However, thanks 
to the D-term instability, // <C 1 might be achieved if Eh is much less than ~ 0.1. Since the 
flat direction completely decays before /-balls are formed for small enough eh, there must 
be an intermediate value of Eh that leads to fu ~ 0.02. 

On the other hand, the new-type Q-balls |7] themselves can be the dark matter in the 
gauge-mediated SUSY breaking models. Since the new-type Q-balls are formed where the 
gravity-mediated potential overwhelms the gauge-mediated one, the above-mentioned argu- 
ments can be applied. If the energy per unit charge of the new-type Q-ball is smaller than 



2 Here we assume that the charge emitted from /-balls due to the Z?-term instability dominate over that 
evaporated from the surface of Q-balls through the interaction with thermal plasma. In fact, the emitted 
charge can be as large as that of Q-balls due to the Z?-term instability. 




(41) 




(42) 
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FIG. 6: Evolutions of the axial ratio of the orbit inside a /-ball with the initial valuers/ = 0.5. 
Due to the D-term instability, ej tends to approach 1 (solid line), while it remains constant in the 
single- field parametrization (dotted line). 



TABLE I: Numerical results of the final axial ratio and /q. For ej < 0.2, the /-ball divides 
into mutliple Q-balls, and the final axial ratio and /q of the largest one are shown in parenthesis 
below. For ej = 0.9, the conversion of the /-ball to a Q-ball cannot be observed in the numerical 
calculation. 
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the nucleon mass, they are stable and can be the dark matter. The relation between the 
dark matter and the baryon abundance is thus given by 

n DM = -\^n B , (43) 

1 - Jb m n 

where we take m<j> ~ m 3 / 2 . This gives the constraint on 



~ 1.0, 0.98, 0.83, for m 3/2 = 0.01, 0.1, lGeV, (44) 

which can be explained by taking e.g., Eh ^ 0.4 if fj ~ 1 and Eh — Si- 
lt has been thought that almost all the baryon number is absorbed in Q-balls. That is 
why it was necessary to estimate how much the baryon charge evaporates from the surface of 
Q-balls through the interaction with thermal plasma. However, in fact, some finite fraction 
is emitted from /-balls due to the D-term instability, and the fraction is determined only 
by the initial axial ratio, independent of the reheating temperature and so on. Thus the 
D-term instability may revive the appealing solution for the coincidence problem between 
the dark matter and baryon abundance, which was otherwise difficult to be achieved. 
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FIG. 7: Spatial distributions of energy and charge of the /-ball with the initial axial ratio: ej = 0.5. 
The /-ball becomes more Q-ball-like by emitting both the energy and the charge. Note that the 
energy is more efficiently released than the charge, which makes the effective axial ratio increase. 
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FIG. 8: Power spectra of the constituent field $i. The model parameters are same as before except 
K\ = K2 = —0.01. The upper panel represents the power spectrum of the initial /-ball solution 
with ei = 0.5, while the D-term instability begins to grow in the lower panel. The instability bands 
coincide with those obtained by the analysis of linearized fluctuations (see Fig. 0). The instability 
bands corresponding to k/m$ = 2,4, 6, . . . are enhanced due to"rescattering" processes. 



IV. FATE OF FLAT DIRECTIONS 

In the previous section we have considered the case that Q-balls are formed. What if 
Q-balls are not formed? If the coefficient of the one-loop correction K is positive, Q-balls 
are not formed. Also even if K is negative, the strong tachyonic instabilities will prevent 
Q-balls from being formed, when the initial axial ratio is small enough. It should be noted 
that Q-balls remain stable under the effect of the D-term instability since the instability 
disappears when the orbit becomes circular. In other words, unless Q-balls are formed, the 
D-term instability will never cease to grow. Therefore, if K < or K > with eh <C 1, 
the flat direction experiences the D-term instability until it deforms into a disordered state, 
which will decay and/or interact with other particles and attain the thermal equilibrium 
sooner or later. The fate of the flat directions is thus very simple: Q-ball or nothing. 
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From the above discussion, the gravity-mediation type Q-balls are formed when the initial 
axial ratio of the orbit is large enough eh > 0.1. Then what about the gauge-mediation 
type Q-balls? Since the strength of the A-term is proportional to the gravitino mass, the 
initial axial ratio tends to be much smaller than 1 if the AD field start to oscillate in 
the gauge- mediated potential Since the homogeneous motion in the gauge- mediated 



potential cannot be solved analytically, it is necessary to perform numerical calculations to 
investigate the D-term instabilities. Due to this additional instability, the Q-ball formation 
will be considerably delayed or prevented. The detailed discussion on this subject will be 
presented elsewhere |20^ . 



V. DISCUSSIONS 



In this paper we have investigated the dynamics of the flat directions taking account of 
the D-term potential. To this end, we adopted the multi-field parametrization of the flat 
directions, which enables us to deal with the relative fluctuations between the constituent 
fields. It is found that there exist instabilities due to the D-term potential and the nature of 
these instabilities depends on the axial ratio e of the oscillating flat direction. If e is close to 
1, the D-term instability is similar to narrow resonance described by the Mathieu equation 
with 5 < 1. On the other hand, the instabilities show the tachyonic nature for e < 1. 
Due to this tachyonic instability, the cosmological evolution of the flat direction are always 
under the influence of D-term instability. In particular we have shown that the existence of 
the instabilities is crucial to the formation of Q-balls. Since the D-term instability violates 
the conservation of the adiabatic charge, /-balls are no longer quasi-stable, and therefore 
the transition from /-balls to Q-balls occurs very efficiently It is important to note that 
during this process some fraction of the charge of the Q-ball is emitted, which may revive the 
relation between the baryon number and dark matter density of the universe. Furthermore, 
the D-term instability drastically changes the evolution of the flat directions. Unless Q-balls 
are formed, the flat direction completely decays relatively soon after it begins to oscillate. We 
believe that the existence of such instabilities is generic to the system of scalar fields with the 
D-term potential. For instance, in the reheating stage of the D-term inflation models |2l| . 
the scalar fields are oscillating in the similar potential. Thus the reheating processes of 
inflation models of this type might be affected by the D-term instability. However, further 
investigation is clearly necessary and is beyond the scope of this paper. 
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